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The generalized Gibbs ensemble introduced for describing few body correlations in exactly solvable 
systems following a quantum quench is related to the nonergodic way in which operators sample, 
in the limit of infinite time after the quench, the quantum correlations present in the initial state. 
The nonergodicity of the correlations is thus shown analytically to imply the equivalence with the 
generalized Gibbs ensemble for quantum Ising and XX spin chains as well as for the Luttinger model 
the thermodynamic limit, and for a broad class of initial states and correlation functions of both 
local and nonlocal operators. 
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I. INTRODUCTION 

Following a number of groundbreaking experiments 
with ultracold atomic systems [IH3] , the problem of ther- 
malization of exactly solvable quantum many-body sys- 
tems has attracted a great deal of attention [JHE]. This 
is because it relates to fairly fundamental questions, such 
as the emergence of thermodynamics in isolated systems 
prepared in initial states that are not eigenstates of the 
Hamiltonian (i.e. systems undergoing a so-called 'quan- 
tum quench'). The latter subject has deep ramifica- 
tions, both in condensed matter physics [7J [5] and cos- 
mology 0E]. Moreover, this problem is also relevant to 
the ongoing efforts to build 'quantum emulators', that 
is, tunable quantum systems capable of accurately sim- 
ulating the mathematical models of many body physics. 
In this regard, the problem of thermalization impacts on 
questions such like how much memory will the emulator 
retain of its initial conditions and whether the standard 
Gibbsian ensembles can be used to predict the outcome 
of the simulation [Jj [8] . 

Interestingly, it was first conjectured by Rigol and 
coworkers [3] that the steady state of simple few-body 
observables of integrable systems following a quantum 
quench can be described by a generalized Gibbs ensem- 
ble (GGE). The density matrix of the GGE is obtained 
as the less biased guess [TU] of the steady state given the 
constraints on the dynamics stemming from the existing 
set of non-trivial integrals of motion. Surprisingly, it was 
found that, in order to reproduce few body observables, 
only a subset of the order of L (where L is the system 
size) of simple integrals of motion is needed [4]. 

Concerning the general applicability of the GGE, there 
has been also some debate [IT] about the importance for 
thermalization of the locality of the operators in the ba- 
sis of eigenmodes of the system. For the quantum Ising 
chain, it was recently shown analytically that correla- 



tion functions of nonlocal operators also thermalize to 
the GGE [12]. Similar results had been found earlier for 
the Luttinger [5] and sine-Gordon models [T5] . 

However, the reason why the GGE has been so suc- 
cessful in explaining the steady-state correlations in very 
diverse models has remained rather obscure. For local 
operators in certain integrable held theories, solid argu- 
ments in favor of the validity of the GGE have been put 
forward by Fioretto and Mussardo [13] • Furthermore, 
Cassidy and coworkers [TS] recently introduced a gen- 
eralization of the eigenstate thermalization hypothesis 
(ETH) [ini H7J for integrable systems. Previously, the 
ETH has been successfully used to understand thermal- 
ization in non-integrable systems [16] . 

In this work, we describe a general method to demon- 
strate the applicability of the GGE in exactly solvable 
models for a general class of initial states. Our method 
does not require the explicit evaluation of correlation 
functions at asymptotically long times after the quench. 
Instead, it suffices to show that the asymptotic correla- 
tion functions of (either local or nonlocal) operators de- 
pend only on the expectation values of quasi-particle oc- 
cupation operators in the initial state. This property, to- 
gether with certain properties of the class of initial states 
considered in this work [531 H3]> allows to demonstrate 
that each eigenmode of the system is subject to a (mode) 
dependent effective temperature, the latter being noth- 
ing but a restatement of the GGE conjecture. This new 
point of view on the GGE also explains some of the less 
well understood aspects of the conjecture that have been 
briefly mentioned above. As a matter of fact, it explains 
why the only set of integrals of motion that are needed 
to construct the GGE correspond to the quasi-particle 
occupation operators. For the class of exactly solvable 
models discussed below, the latter are a minimal set of 
L integrals of motion that entirely determine the asypm- 
totic correlations. The fact that the asymptotic corre- 
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lation functions depend only on the expectation value 
of these non-trivial integrals of motion means that the 
system remembers much more information about its ini- 
tial conditions than it is the case in systems exhibiting 
thermalization to a standard Gibbs ensemble. For the 
latter, only the expectation value of the energy and the 
Hamiltonian suffice to determine the effective tempera- 
ture and chemical potential of the standard Gibbs en- 
semble describing thermal equilibrium of systems in the 
thermodynamic limit. The lack of relaxation of correla- 
tion functions to thermal equilibrium found in this work 
bears a strong resemblance with the nonergodicity of the 
magnetization in the XY model discussed several decades 
ago by McCoy [18] and Mazur [19] (see [8] for a recent 
review of this result). Thus, we shall call this behavior 
of the asymptotic correlations 'nonergodic'. 

Our goal in this article will be to illustrate our method 
to demonstrate the applicability of the GGE by applying 
it to several models that have been previously analyzed 
either analytically [51 ITS] [20] or numerically [31 [15]: The 
quantum Ising chain (cf. Sec. |TT[ ) , the Luttinger model 
(cf. Sec. Ill), and the lattice gas of hard-core bosons in 



one dimension or quantum XX spin chain (cf. Sec. IV). 



In the former two cases, we consider a quench from an ini- 
tial state that does not break the (lattice) translational 
invariance and it is therefore conceptually simpler. In 
Sec. |IV[ we turn to a the more involved case in which the 
initial state is not translationally invariant but conserve 
the particle number. In our analysis, we shall consider 
correlations of both local and nonlocal operators, demon- 
strating, for a broad class of initial states, that in both 
cases thermalization to the GGE takes place. We also 
consider a more general type of initial states that those 
that have been investigated in the past [5J HU [20] ■ We 
have relegated to the appendices the discussion of some 
of the most technical aspects of this work. 



II. THE QUANTUM ISING CHAIN 

Let us begin with an precise statement of the prob- 
lem that we intend to address. A quantum quench 
refers to the situation in which a system is prepared 
at t — in an initial state (denoted po below) that is 
not an eigenstate of the Hamiltonian H. Furthermore, 
we shall assume that, following the quench, the system 
reached some kind of equilibrated state where observables 
and correlation functions acquire (time-averaged) values 
about which they exhibit small temporal fluctuations (cf. 
Fig. [IJ. A necessary condition for this equilibration to 
occur is that if we expand po in the basis of eigenstates 
of H, 



Po 



(i) 



the coefficients C mn are sufficiently nonsparse on the ba- 
sis of eigenstates of H such that (for a sufficiently large 
system) unitary evolution can lead to a steady state as 



a result of dephasing between the contributions of many 
different eigenstates to the expectation value of observ- 
ables and correlation functions. The general conditions 
for equilibration have been discussed in Ref. [2"TJ . 

We next begin our investigation of quantum quenches 
in exactly solvable models by considering the quantum 
Ising chain. For this system, the Hamiltonian that de- 
scribes the time-evolution of the system following the ini- 
tial state preparation takes the form: 



(2) 



where erf and erf are the Pauli matrices at site j and 
J and g are the model parameters. As reviewed in Ap- 
pendix|Aj the Hamiltonian in Eq. Q can be diagonalized 
by means of a non-local transformation due to Jordan 
and Wigncr, which uncovers the fact that its elementary 
excitations are indeed free fermions described by 



H = ]>>#) 



(3) 



where e g (k) = 2 J-y/ 1 + g 2 — 2g cos k is the fermion dis- 
persion (|fc| < 7r). The operators 'y(k) are the eigenmodes 
of the system for the set of parameters (J, g). They evolve 
according to the law: j(k,t) — e im l h ^[] t ^ e r im l h — 
e -««s( fe )*/fi/y(^). However, the actual observables of the 
system correspond to the Pauli matrices, of and erf. For 
example (cf. appendix |A[), 



erf = 1-2/1/ 



where 



jJ3> 



— ie g (k)tfh 



(4) 



7 (fc) 



+<.(x)e+"' (fe)t/ V(fc) 

= sin k I (cos k - 
] cos{9 g {k)/2)/VZ and 
such that e.g. 



(5) 



9 

v k {x) 
\u k (x)\ 2 



and 



where tan0 s (fc) 

Wfe(x) = 

-ie ikx sin(0 s (fc)/2)/V% 
\v k (x)\ 2 = L-\ 

The class of initial states with which we shall be con- 
cerned in what follows is described by a density matrix 
Pa = Zq 1 e~ H °/ T , where Zq is a normalization constant 
and the operator Hq is a quadratic form of the eigenmode 
operators j(k) and 7^(fc) (cf. Eq. [6]). Since Hq must be 
hermitian, it can be interpreted as the Hamiltonian of 
the system at time t < 0, and the parameter T as the 
absolute temperature of an energy reservoir with which 
the system was in contact (the pure state case is obtained 
by taking T — > 0). The contact with the reservoir is re- 
moved at t — as the Hamiltonian is suddenly changed 
from H to H and the system allowed to evolve unitarily. 
This defines the kind of quantum quench that has been 
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analyzed in most cases so far [3HS1 QU HI- Thus, the 
most general form for the initial Hamiltonian Hq reads: 



H = [eo(k)8k,k> + V (k, k')] 7 f (fc)7(fc') 

k,k' 

+ [AS(fc, k)j(k)j(k') + A (k, fc / ) 7 t (*')7 t (*)] ■ (6) 

k,k' 

The term proportional to Vo(k^k') in Eq. ^ can be in- 
terpreted as a scattering potential that is switched off at 
t = 0. The presence of the scattering potential in Hq 
means that, in general, the initial state, po, breaks the 
translational invariance of the lattice. Furthermore, the 
last two terms in Eq. Q imply that that the number of 
fermion quasi-particles is not well defined in the initial 
state because [po,N] ^ 0, where N = X)&7 1 is 
the quasi-particle number operator. 

The states introduced above have two important prop- 
erties: i) Correlations of products of an arbitrary (even) 
number of Fermi operators like j(k), J^(k), or and 
fj, can be expressed in terms of products of correla- 
tion functions of bilinear operators like (k)j(k')) = 
Tr po7 (k)'y(k'), (7(fc)7(fc')), etc. This result is known as 
Wick's (or more precisely, Bloch-de Dominicis' [22]) the- 
orem and it is needed to show that the correlation func- 
tions of nonlocal operators can be obtained from those 
of local operators (see below and Appendix |B| ; ii) for 
any partition of the cigemodes into two disjoint subsets 
A (called "system" in what follows) and B (called "envi- 
ronment"), the reduced density matrix obtained by trac- 
ing out the environment B can be also written as the 
exponential of a quadratic form of the Fermi operators 
7(fc) and 7 t (fc) [23J 123 • In particular, if the "system" A 
consists of a single eigenmode (and therefore B contains 
the remaining L — 1 modes), the reduced density matrix 



p(k) = Tik'^k Pa = Z 1 (k)e 



-X(k)I(k) 



(7) 



where the symbol Tv k i^ k stands for the partial trace over 
all modes but k and I(k) = ^(k)^(k) is the quasi-particle 
occupation operator. This result applies to systems with 
eigenmodes obeying both Fermi and Bose statistics 
24'. It is worth noting that the GGE density matrix is 
constructed as the direct product of these single-mode 
reduced density matrices: 



Pgge 



)p(k), 



(8) 



which is nothing but the mathematical statement that 
each mode is subject to a mode-dependent effective tem- 
perature T(k) = X(k)/e(k). In fact, our method to prove 
the applicability of the GGE will rely on this interpreta- 
tion of the GGE. 

In order to make contact with previous studies of the 
quantum Ising chain [11] [T21 [20] , we shall analyze the fol- 
lowing the case of an initial state that respects the lattice 
translation symmetry. This requires that Vo(k,k') = 



and A (k,k') — ^h (k)S k+k > ^- For specific choices 
of eo(fc) and ho(k) such that eo(k)/ho(k) = tan^(fc), 
with cj>(k) = 6 g {k) — 6 go (k), H Q would correspond to 
the Hamiltonian of the quatum Ising chain at a dif- 
ferent value of the parameter g — go- At such point, 
e so(^) = \/ e o(k) + h-o{k), is the dispersion of the quasi- 
particles which are described by a different set of eigen- 
modes related to 'y(k) and 7' (k) by a canonical transfor- 
mation parametrized by the angle <fi(k) [5D]. Thus, for 
such a choice we can speak of a quench in the parameter 
g. However, for arbitrary eg(k) and ho(k), H does not 
map to a quantum Ising chain Hamiltonian, and there- 
fore, our choice of the initial state of the quench, albeit 
translationally invariant, is more general than previous 
choices [TTJ [TJJ [2U] , which focused on quenching the pa- 
rameter g only. 

We next turn to the correlation functions of the model 
following the quench. We begin with the discussion of 
the correlation function for a local operator such like the 
Fermi field: 

cf){x i -x j ,t) = {f}{t)f j {t)) (9) 
= Y, [«*(a:i)«fc'(^)e < [ e « ( * ) - e ' (fc ' ) ]* / ' l <7 t (*)7(A , )> 

k,k> 

+« fc (a; l )^(x J ) e - l h«- £ 9 ( fe ')] t / ;i (7(fc)7t(A ; '))" 
+ ^[<(x l )^(x J )e 4 ^W + ^( fe ')]*/ ri ( 7 t(fc)7(fc')) 

k,k' 

+ Wfc (x l )u fc '(^)e- l ^ (fc » +e '(' £ '^/ ri (7(fc)7 t (fc')>] • (10) 

Thus, at any finite t and for an arbitrary initial state, 
the above correlation function is fully determined by the 
eigenmode correlations in the initial state Go(k,k') = 
(7 t (fc) 7 (fc')) and F (k,k') = (j(k)-f(k')) . However, the 
invariance of the initial state with respect to lattice trans- 
lations greatly simplifies the above expression implying 
that G (k, k') = S kM N (k) and F {k, k') = F (k)5 k+ k>,o- 
Hence, 



f {Xi Xj , t) 



k 



[ul(xi)uk(xj)N (k) 



+v k (x i )v* k (x j )(l + N (k))\ 
[ut(x i )v k ( Xj )F (k)e 2i ^/ H (11) 



+u k (x l )v* k {x J )F*{k)e- 2 ^ t / h 



(12) 



We shall next consider the limit t — > +00 of the 
above expression after taking the thermodynamic 
limit where L — > +00. Note that u* k {xi)v k {xj) — 
-iL- 1 e ik(x >- x ^sin(6 g (k)/2)cos(6 g (k)/2) and F (k), 
which is itself a function of eo(/c) and ho(k) are assumed 
to be well behaved, smooth functions of k. Therefore, 
it follows, by virtue of the Riemann-Lebesgue lemma, 
that the second term in the right hand side of Eq. ( 12 ) 
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vanishes in the t 



Df}(xi — Xj) = lim Cfffe — Xj, t) 
= ^2 [ul(xi)u k {xj)N (k) 



-00 limit. Thus, 



(2), 



(13) 



+v k (x i )v* k (x j )(l-N (k))], (14) 

where the thermodynamic limit is implicitly understood. 
Note that the above result, Eq. ( fi"4| , means that this 
correlation function depends only on the expectation val- 
ues of the L integrals of motion I(k) = y{k)j(k). In- 
(2) 

deed, D\J(xi — Xj) is a (weighted) sum of the expec- 
tation values, N (k) — (I(k)) — Trp I(k). Hence, for 
each term of the sum over k, we can use the second of 
the properties of the class of states po described above, 
namely, we can trace out all the modes k' 7^ k and write 
N (k) = Tr p(k)I(k), where p(k) is given in Eq. 0, 
with X(k) = In [(N (k) - 1)/N (k)}. This result obtained 
via a partial trace amounts to the statement that each 
eigenmode is subject to a mode-dependent effective tem- 
perature, which is equivalent to conjecturing that the 
asymptotic state is described by the GGE density matrix, 
Eq. pj . This result also implies that the Cff(xi — Xj , t) 
will not relax its thermal equilibrium value, a behavior 
that we call 'nonergodic' [HI \T§\ . 

Similar results can be obtained for the asymptotic limit 
of other correlation functions of local operators such as 



.4 



/] + fj and Bj — f] — /,-. We merely state here 



J J 3 

the results 

>(2). 



D AA(xi ~Xj) = t Hrn^ C%X(xi - xj,t) 

D%{xi - Xj) = ^ Urrx^ C%{xi - xj,t) (17) 
= lim (BMBjQ)) = -6ij (18) 



(15) 
(16) 



Dab( x i ~ x j) — . 1™ ^ief^ — x ji~k) 

t—$-\-OQ 

= lim {Mt)B 3 {t)) 

t— >-\-OD 



(19) 



+ e<*(»«-*i> e <0 «W(l-iVo(ife))]. 



(20) 



Again we find that the asymptotic correlation functions 
are nonergodic, as they only depend on No(k). 

Using the above results, we are now in a position to 
discuss the correlations of a nonlocal operator such as 
erf . Nonlocal means that this operator does not reduce to 
a simple linear combination of the eigenmode operators 
7(fc) and 7^(fc). Indeed (cf. Appendix |A| , 



(/, 



Oil* 1 



(21) 



that is, erf involves an infinite product of local operators 
(in this case fi and fj). As it is discussed in the Appen- 
dices [B] and [C] the two-point correlation function of erf 
can be expressed, by means of Wick's theorem, in terms 
of a finite product of (equal time) correlation functions 
of the local operators A* and Bi. The existence of the 
t — > +00 limit of those correlation functions (cf. Eqs. 16 



18 20 ) suffices to ensure the existence of the asymptotic 



correlation function (cf. Appendix O): 



D^(x i -x j )= firn (a* - x 3 , t) (22) 

t— >- + oo 

= MK(^i. (23) 

t— f +00 J 



In this limit (a thermodynamically large system is im- 
plicitly assumed) , the above correlation function reduces 
to a finite Toeplitz determinant (see Appendix [C]) which 
depends on Dy B (x k —x{) (cf. Eq. p0| ) with i < k,l < j. 

Thus, we conclude that just as for the local correla- 
tions discussed above, the nonlocal correlations are also 
nonergodic and are given by the GGE, which assumes a 
mode-dependent effective temperature. As a corollary, it 
also follows that only the L occupation numbers No (k) = 
(7 WtW) are needed to determine the asymptotic cor- 
relations of both local and nonlocal operators. Other 
integrals of motion different from I(k) = (k)"f(k), such 
like e.g. the products I(ki)I(k2), etc. do not play a role 
in determining the asymptotic correlations and in the 
GGE. The set of L occupation numbers, N (k), amounts 
to much less information than the full initial-state corre- 
lations, which are determined by both iVo(fc) and Fo(k) 
(3L real numbers, in total). However, these L occupa- 
tion numbers amount by far to much more information 
that the expectation value of the energy (H) = Tr po H 
and the particle number (N) , which determine the effec- 
tive temperature and chemical potential in the case of 
thermalization to the standard (grand canonical) Gibbs 
ensemble. 

In this section we have focused on the quantum Ising 
model, which exhibits fermionic quasi-particles. How- 
ever, this is not a limitation to our methods, as shown in 
the following section, where we deal with the Luttinger 
model exhibiting bosonic quasi-particles. We have also 
required that the initial state respect lattice translational 



j<i 



invariance. As we show in Sec. IV this is again not a se- 
rious limitation to demonstrate the applicability of the 
GGE. 



III. QUENCH IN THE LUTTINGER MODEL 

Let us next consider a quantum quench in the Lut- 
tinger model (LM) [SJ US] , which is a model exhibiting 
bosonic quasi-particles. The initial state is assumed to 
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be of the form po oc e Hq / t , where 



ifld(fc)[6 t (fc)6 t (-fe)+6(fc)6(-fc))]}, (24) 



vo(k) and go (A;) being regular functions at k = 0. 
The operators 6(fc) and 6^(fc) obey Bose statistics: 
[b(k), b'(k')] = 6k,k'j commuting otherwise. They are 
eigenmodes of the Hamiltonian 



H = J2 hv ( k )\ k \ b f (k)b(k), 



(25) 



which dictates the evolution of the system for t > and 
which we assume to describe an interacting version of 
the LM. Differently from the XX chain studied in the 
previous section, the eigenmodes of the LM are bosonic. 
In the initial state po, Eq. (24), the number of bosonic 



modes is not well defined since 



[fT , $>*(*)&(*)] ^0. 
Mo 



(26) 



However, po does commute with the momentum opera- 
tor P = ktf (k)b(k) , which implies that that po is a 
translationally invariant state. 

We shall assume below that the 'fundamental' fermions 
of the model [SJ [2^] also diagonalize Ho. This amounts 
to assuming that Hq describes a non-interacting version 
of the LM [5]. Therefore, we can regard this situation as 
a quench from the non-interacting to the interacting LM, 
with interactions V = H — Hq being suddenly turned on 
at t = as the contact with a bath at absolute temper- 
ature T is also removed [5]. This allows us to determine 
the relation of the eigenmodes to the observables of the 
system. 

Physical operators in the LM can be expressed in terms 
of exponentials or derivatives of (chiral) boson fields de- 
fined as follows: 



<j> a {x,t) 



_ iHt/h 



4> a {x)e 



iHt/h 



00a (*) + ~l~ N <* 



■* a {x,t) + & a (x,t), (27) 



where L is the system size and N a is the number of 
fermions of chirality a — r, I ([N a ,<f>op] — iS a p). Be- 
low, we shall work in the sector where N r = Ni = 
which also contains the ground state of H, namely |0) 
(i.e. b(k)\0) = for all k). Furthermore, in terms of the 
eigenmodes [51 |2T)]. 



$ a (x,t) = (fr) 2 elSakx [cosh0(fc) e - TO(fe)|fe|t 

x b(s a k) - e m{kmt sinhO(k) b\-s a k)] , (28) 



with s r — — si = 1 and tanh0(£;) = g (k)/vo(k). Using 
these chiral fields, the density (or 'current') operator for 
each fermion chirality reads 

J a (x,t) = : ^ a (x,t)ip a (x,t) : = ^-d x (/) a (x,t), (29) 

Zir 

: . . . : meaning normal order with respect to the ground 
state of H [53]. Note that the J a (x,t) are local in the 
eigenmodes, b(k) and tf(k). On the other hand, the 'fun- 
damental' fermion fields 29 , 



4>a(x,t) 



is a cj> a (x,t) 



(30) 



are nonlocal ('vertex') operators. Using Wick's theorem, 
we can recast any fermion correlation function in terms 
of two body correlators of the local fields (j> a (x, t) because 
the cumulant expansion to second order is exact for states 
like pq. Mathematically, 



a iA a (xi,—,x n ,t)\ 



e -|<A=(xi,...,x„,t)) 



where Q^Pi = 0) 
A a (xx, ■ 



, X. 



i,t) = / ]Pi(j>a{Xi,t). 
i=l 



(31) 



(32) 



Eq. (31) can be proven by expanding in series the ex- 



ponential in the left hand side in a Taylor series and 
applying Wick's theorem to all the terms, which involve 
powers of A(xi, . . . ,x n ). Resuming the resulting series, 
the right hand side of pTl) is obtained. 



From the previous discussion, it can be seen that in 
order to compute the equal time correlations of the LM 
it is sufficient to consider 



cgW) = (M^t)H0,t)) = D%(x)+F£>(x,t), (33) 
where 



p(2), 



7T 

\k~\L 

ik 



(cosh26»(fc) + sgn(fc)) 



x [e lkx (1 + N (k)) + e- tkx N (k)] , (34) 

is the contribution of the diagonal correlations No(k) = 
(&t(jfc)6(fc)). However, 



(cosh26>(fc) + sgn(fc)) 



Jkx-2iv(k)\k\tp Q ^ + e -ikx+iv{k)\k\t F *(j^ ^ ^ 



where Fo(k) — (b(k)b(-k)) , stems from the 'anomalous' 
(i.e. 'superfluid') correlations. Note that the transla- 
tional invariance of the initial state implies that No(k) 
and Fo(k) are the only non- vanishing two-point correla- 
tions of the eigenmodes in the initial state. Whereas the 
contribution of the diagonal correlations N (k) is time 
independent, the contribution of the anomalous terms 



G 



depends on time. It may be expected that, because 
of dephasing between the different Fourier components 
(i.e. the Riemann-Lebesgue lemma), in the thermody- 

( 2^ 

namic limit (x,t) vanishes as t — > +00. However, 
the t — > +00 limit of this function limit must be handled 



operators fi , fj and Fourier expand the latter in terms 
of /(fc) and P(k) (cf. Appendix |A| . Hence, 



/r = 5>(fc)/t(*o/(fc), 



(39) 



with care because the l/|fc| in Eq. (351 yields terms di- 
verging logarithmically as t — > +00 5 . Fortunately, as 



we have described above (cf. Eqs 29 31 1, only the deriva- 

(2) — — 

tives or exponentials of C , (x, t) appear in the physical 
correlation functions of the LM model. For example, us- 



ing (30), the two-point correlation function of the right z~ x er H °l T where 



with e(fc) = — Jcosfcao — h (0 < fc < it). Thus, the 
eigenmodes of the system are described by the Fermi op- 
erators /(fc) and /^(fc), which evolve in time according 
to /(fc, t) = e iHt / h f{k)e~ lHt / h = e- te< - k ^/ h f{k), etc. 
The initial state is given by the density matrix po = 



moving Fermi fields reads: 



(V40M)Vv(O,t))=Aexp cf\ Xl t)-C^(0,t) . (36) 



.(2), 



Since (for L -» +00 and T = 0) 
J^(x,t)-i^(0,t)~log 



(2vt) 2 - x 2 



(2vt) 2 



(37) 



where v = v(k = 0), the time-dependent logarithmic con- 
tributions disappear (for finite x) as t — > +00 [5]. Thcrc- 
fore, we can safely ignore the contribution of ' [x, t) 
in the t —> +00 limit. This means that all correla- 
tions are asymptotically determined by D<p r (x), which 
only depends on N (k) — (tf(k)b(k)), i.e. it is non- 
ergodic. Furthermore, for each term of the sum in 
Eq. (34), we can trace out all the fc' 7^ fc and write 
N (k)^ Tr [p(k)tf{k)b(k)] = Tr p GGE tf(k)b(k). Since 
p(k) = Tr k ^ k , p = Z- X [k) e -Hk)bUk)b(k) ^ we arrive 
at the same result as if we had used the GGE density 
matrix pgge = ®fc/°(^)- Thus the equivalence with the 
GGE is established for the simple correlation funcitions 
involving the bose field <p a (x) in the LM. 

Finally, it is interesting to note that translationally 
invariance requires that eigenmode correlations are bi- 
partite, that is, each mode fc is correlated only with the 
eigenmode at — fc (cf. Eq. (24)). Thus, an alternative 
way of obtaining the results of this section and those of 
Sec. [n] is to compute the reduced density p(k) as a partial 
trace for a partition of the eigenmodes into fc > and 
fc < 0. Thus, we can regard the effective temperature 
T(fc) = \(k)/v(k)\k\ for e.g. the modes with fc > as the 
result of their correlations with the fc < modes (and 
viceversa) [2"3"] . 



IV. THE XX CHAIN 



The Hamiltonian of the XX chain reads 

i=l i=l 

in terms of the Pauli matrices af , of . We shall assume 
an open chain like in Ref. [4]. In order to diagonalize 
the Hamiltonian, we first carry a Jordan- Wigner trans- 
formation to express the Pauli matrices in terms of Fermi 



H = J2 [eo(*)*M' + Vo(k, fc')] f\k)f{k'). (40) 

k,k< 

In order to make contact with the numerical studies of 
Ref. [J] , we have assumed that the initial state commutes 
with N in Eq. (40). Therefore, the anomalous terms 
(such like those cx Ao(fc, fc') in Eq. [6]) are absent in this 
case. However, the presence of the scattering potential 
Vb(fc, fc') implies that the initial state po breaks the lattice 
translational invariance. 

We next turn to the analysis of correlation functions. 
We first consider the equal time correlation of a local 
operator like 0(xi) = J^kV'kixi) /(fc), where we shall 
require that (the square of) fk{xi) is normalized to the 
system size (i.e. |(/?/c(xi)| ~ 0(L _1 / 2 )). This means that 
the quasi-particles of the system (described by the eigen- 
modes /(fc) and P (fc) of the Hamiltonian H) can occupy 
extended orbitals after the quench and are not localized. 
In other words, in the thermodynamic limit the quasi- 
particle spectrum of H is assumed to be described by a 
continuum of extended (i.e spatially delocalized) levels 
with no macroscopic degeneracies. In principle, violation 
of this requirement may prevent the system from reach- 
ing equilibration as contributions from localized states 
will lead to oscillatory behavior at long times and the 
absence of decoherence. With this caveat, let us con- 
sider: 

C [ o\x u x 3 ,t) = (0^( Xi ,t)0(xj,t)) 

= ¥>fc(*i)M*i) Go(K fc') e -^«-^')]*A (4i) 

k,k' 

which depends on the eigenmode correlations Go(fc, fc') = 
(P (k)f(k')). The latter are L 2 real numbers containing 
the full information about the initial state [23 [23] . With 
the above assumptions and in the thermodynamic limit, 
we find that (see discussion below) 



(2)^. _ lim C (2)( 

Y / 'fil(^Wk(x J )No(k), (42) 



Dy(xi,Xj) = ^ lirn^Q, (x i} Xj,t) 



where N (k) = (/(fc)) with /(fc) = 7^ (k)'y(k) is the quasi- 
particle occupation in the initial state. Eq. (|42|) means 

that D (xi,Xj) is nonergodic [HI HE] as it is entirely 



7 




1000 2000 3000 4000 5000 



u/h 



FIG. 1: (color online) Time evolution of the correlation func- 
tion of the local operator O(xi) = f, for i — L/2 and L/2 + 1, 
(J$(t)fi+x{t)), for L = 800 and L = 1600, with TV = 40 and 
TV = 80 . The horizontal line corresponds to the time average 
for the largest system size (L = 1600 and TV = 80). In the 
initial state TV = 40 (TV = 80) hardcore bosons that are con- 
fined by a harmonic potential of the form Vo(xi — xm) 2 , with 
Vb = 10~ 3 J/y/N and x M = L/2. 



determined by L real numbers, the quasi-particle occu- 
pations N (k) — Tr \p(k)I(k)]), where p(k) has been de- 
fined above (cf. Eq. [7|. Hence, we can again perform the 
partial trace in each of the terms of the sum ( 42 ) , and 



conclude that each eigenmode is subject to a fc-dependent 
effective temperature, as expected from the GGE. This 
result also implies that the asymptotic correlation func- 
tions are determined by much less information than the 
one contained in the initial state (i.e. O(L) vs. 0(L 2 ) 
real numbers). Yet, this is much more information than 
the one needed to characterize the asymptotic state of 
standard thermal equilibrium. 

In order to demonstrate Eq. ( 42 ) , we display in Fig. [I] 
the results of the numerical evaluation of the time evo- 
lution of Cf)(xi,Xi+i,t) (i — L/2) using Eq. ( [41] ) for 
O(xi) — fi. We consider an initial state po for which 
Vo(k, k') in Eq. (40) is a harmonic potential that confines 
TV hardcore bosons at the center of an open chain of L 
sites. The potential Vp fk, k ') is taken to scale as Vq/vN 
(Vo = 10~ 3 J in Figs. Tj2 and [5J in order to obtain a 



well defined thermodynamic limit of the initial cloud of 
hardcore bosons |T5] [25] . This potential is switched off 
at t = and the bosons are allowed to expand [4]. The 
vertical line in this figure corresponds to the time average 
(for L = 1600, the average for L = 800 is not shown but 
it is very close to it). The average is given by Eq. ( [43] ) 
evaluated at finite L. It can be seen from Fig. [T] that 
for both L = 800 and L = 1600, after a short transient, 
the correlation function exhibits roughly equilibrationa 
and its time fluctuations about the average become fairly 
small. In can be also seen that, as L increases from 800 to 
1600 (while keeping N/L constant), the size of the time 
fluctuations decreases suggesting that for L — > +oo they 
will vanish. Therefore, in the thermodynamic limit the 
asymptotic correlations are given by the quasi-particle 



5^ 0.04 

-a 

5 0.03 

— 0.02 

II 

ST o.oi 

O 0.00 



Jill 



Vfh = 0.05 

fc'/jr = 0.10 



-k'/x = 0.20 



, miMlllllii . 



k'/ir = 0.40 



-fc'/jr = 0.30 



k/TT 

L = 800 
TV = 40 



0.0 



0.2 



0.4 0.6 

k/TT 



0.8 



1.0 



FIG. 2: (color online) Modulus square of the eigenmode cor- 
relations, |Go(fc, k')\ 2 = j(/ t (fc)/(fc'))| 2 in the initial state of a 
XX chain at T = 0, for a system in a box of size L = 800 con- 
taining TV = 40 hardcore bosons, and fc'/V = 0.05, 0.1, 0.2, 0.3 
and 0.4 (inset). The initial state describes TV = 40 hard- 
core bosons that are confined by a harmonic potential with 
the same parameters as in Fig. [I] Note that the eigenmode 
correlations |Go(fc,fc')| 2 are strongly peaked at k = k! . For 
k > ko with ko/iv ~ 0.35 for TV = 40 (see inset) the correla- 
tions are no longer peaked at k ~ k! . However, they become 
substantially smaller than the typical peak values at k' < ko 
(ko increases with TV, cf. Fig. Hh. 



occupation iVo(fc). In the Appendix [B| we further ex- 
plore the equivalence between the thermodynamic limit 
of time-averaged correlations and their t — > +oo limit 
after taking the thermodynamic limit. 

To understand the behavior displayed in Fig.[l]in phys- 
ical terms, note that, in the thermodynamic limit, the 
sums over k and k in the expression for Cq (Xi,Xj,t) 
become integrals and dephasing between different eigen- 
mode contributions to Eq. (41 1 leads to the decay in time 



of the correlations except for the terms where e(fc) = 
e(fc'). It is worth investigating how this dephasing takes 
place in more detail because, generally speaking, for non- 
translationally invariant states, Go(k,k') is not gener- 
ally speaking a smooth function of k and k' (cf. fig- 
ure [2] and [3J) . Thus, arguments based on the Riemann- 
Lebesgue lemma similar to those employed in sections |n| 
and |III| for translationally invariant states are not appli- 
cable. However, for non-translationally invariant states 
(i.e. Vo(k,k') ^ 0), we numerically find that, as the 
thermodynamic limit is approached (cf. Fig. |3| 



|G (M')| 2 -> Nl(k)5 Kk , + AR (k,k') 



(43) 



where ARo(k,k') decays rapidly for \k — k'\ ^> L^ 1 . 
Eq.(43) must be understood as the statement that typ- 
ical correlations become strongly peaked at k — k' as 
L —> +oo. Thus, setting Go(k,k') ~ TVo(fc)<5fc.fc' in 



Eq. (41 ) becomes an increasingly good approximation at 



large t where decoherence acts most efficiently on the con- 
tributions to the double sum ( 41 ) of quasi-particle levels 



k and k' that are close in energy and correlated (i.e. for 
which Go(k,k') is not negligibly small). 

The claim of Eq. ( 43 ) is illustrated in figures [2] and [3] 
for the same system used to generate Fig. [l] Fig. [2] dis- 
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FIG. 3: (color online) Finite-size scaling of the Mod- 
ulus square of the eigenmode correlations, \Go(k, k')\ 2 = 
\(f (k)f(k')}\ 2 for the same initial state of a XX chain as 
described in the caption of Figs. [T] and [2] for systems of size 
L = 400,800, 1600 and fc'/vr = 0.25. Note that the peak be- 
comes narrower as the system approaches the thermodynamic 
limit where N and L — > +oo while no = N/L = const. 



plays |Go(fc,fc')| 2 as a function of k for several values of 
k', for a system of L = 800 sites and N = J2k N o( k ) = 40 
hardcore bosons at T — 0. It can be seen that \Go(k, k')\ 2 
is strongly peaked at k = k' for k' < ko (kg/ir ~ 0.35 for 
N — 40). However, for k! > ko (see inset) the peak is 
no longer at k = k'. In this case, however, the values of 
\G (k,k')\ 2 become very small compared to typical the 
peak values of \Go(k, k')\ 2 for k < k . The cut-off fc is 
determined by the number of hardcore bosons in initial 
state N oc L (see discussion below and Fig. [4]). As the 
system size L increases while keeping the lattice filling 
?io = N/L constant and scaling the initial trap strength 
as Vq/VN [S3], the peak of Go(k,k r ) at k = k' becomes 
substantially narrower (cf . Fig. [3]) . We also found a simi- 
lar behavior of Go(k, k') when the potential that initially 
acts upon the hardcore bosons was taken to be an ex- 
tended superlattice [55] and an infinite square well box 
of size smaller than L [27] . Although we cannot find 
a rigorous mathematical proof that Eq. (43) holds for 



any non-translationally invariant initial state of the form 
Po ~ e~ H "/ T , with Hq given by (40 1, we expect it to hold 
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FIG. 4: (color online) Log-log plot of the size dependence of 
the cut-off ko for no = N/L = 0.05. The continuous curve is 
just a gide to the eye. Other parameters are the same used 
to generate figures [l] [2] and [3] 



for any physically sensible scattering potential Vo(k,k'). 
The following physical argument can be used to support 
this expectation. Recalling the relationship between the 
eigenmodes of H and the eigenmodes of the initial Hamil- 
tonian, Hq, the two sets of operators are related by the 
canonical transformation, f(k) = ^ a 4>a{k)f(a); 4> a (k) 
are the eigenfunctions of Ho in the basis of eigen orbitals 
of H, that is, 4> a {k) — (k\a) where Hq\cx) = (e Q — fJ-)\&)- 
Hence, 



Go(k, k') = K( k )<P a (k')F(e ai T), 



(44) 



where F(e a ,T) = [e( e ° _A1 ^ T + l] 1 are the thermal oc- 
cupations of the orbitals 4> a {k) in the initial state and p is 
the chemical potential. The latter allows us to fix the av- 
erage number of quasi-particles in the initial state. Next, 
we notice that, as the thermodynamic limit is approached 
(i.e. L — >• +oo while keeping the lattice filling no = N/L 
constant), Eq. (44 1 becomes an infinite sum. Correlations 
are maximized for k = k' because the summands, | (a\k) | 2 
are all positive and correspond to the probability that a 
quasi-particle initially found in the state \a) ends up in 
the state |fc) after the quench. On the other hand, for 
\k — k'\ 3> L , the eigenmodes become largely uncorre- 
lated because the sum in Eq. ( 44 1 involves a large num- 
ber of complex amplitudes (k\a)(a\k') = cf)* a (k)(j) a {k'), 
describing the quantum interference between transitions 
where the quasi-particle ends in one of two orthogonal or- 
bitals, either \k) or in \k'). As L — > +oo, the amplitudes 
interfere destructively and typically average to zero. Fur- 
thermore, at low temperatures, the sum (44) contains a 
cut-off which is roughly given by the state a max for which 
F(eo, max — jj,) <^ 1. Thus, for example, for T = 0, |a max ) 
corresponds to the orbital with the N-th smallest eigen- 
value e a — fi. As k — k! increases, the number of nodes 
of the orbital \k) in the initially (harmonically) confined 
region increases, and the overlaps <j) a {k) = (k\a) dramat- 
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ically decrease in magnitude. This explains the existence 
of the cut-off fco seen in Fig. [2] Consistent with this effect, 
we numerically observe (see Fig. [3]) that k increases as 
the number of hardcore bosons N increases, which makes 
owx bigger. At higher temperatures, the increase in en- 
tropy implies that quasi-particles are spread even more 
over the set of initial orbitals \a), and the quantum in- 
terference effects will be weakened thus decreasing the 
correlations \Go(k, k! ^ k)\ 2 even further. 

To end this section, we shall use the above results to 
sketch the proof that the momentum distribution func- 
tion of the hardcore bosons, S(q,t), for t — > +00 is also 
nonergodic. We first recall that 

S(q,t) = Lj2z iq(Xi - Xj) <^ r (*)4(*)> (45) 



L 



J2e^-^C£Hx llXj ,t), (46) 



(2) 

where Ci a ;(xi,Xj,t) = (af(t)aj(t)). The latter correla- 
tion function involves the operator a x that is nonlocal in 
the eigenmodes f(k),f^(k). However, using Wick's the- 
orem, it can be written in terms of products of two-point 
correlation functions of the local operators A4 = fj + fi 
and Bi = fj — fi. The t — > +00 limit of the product 
exists provided the t — > +00 limit of the two-point corre- 
lation functions involved in the product also exists (see 
the extended discussion in Appendix [B]) . Moreover, as 
t -> +00 



(47) 
(48) 



(2) 

Hence, the t — > +00 limit of C xx \Xi,Xj,t) reduces to a 
Tocplitz determinant involving lim t _j. +00 (yli(t)i3j(<)) = 

D^g(xi,Xj) only (cf. Appendix [c| , where 

Df B { Xi ,Xj) = -Sij +2j2v%(x i )<p k (x j )N (k), (49) 



with tpk(xi) = \/2/{L + 1) sin kxi. This result implies 

that (xi, Xj) is nonergodic. Hence, C x \ (%i , Xj , t — > 
+00) and S(q,t — > +00) are also nonergodic. 

V. DISCUSSION, SUMMARY, AND OUTLOOK 

We have shown that, in exactly solvable models, the 
correlation functions of both local and non-local oper- 
ators, at asymptotically long times, are functions of the 
quasi-particle occupations in the initial state, for a broad 
class of initial states. This means that correlation func- 
tions in these systems retain much more memory of the 
initial conditions than in systems relaxation to thermal 
equilibrium, which is described by a standard Gibbs en- 
semble. This lack of relaxation is similar to the observa- 
tions of McCoy [IE] and Mazur [TH] for the magnetization 



in the XY model. It implies equilibration [ST] but lack 
of ergodicity [19] , as the existence of non trivial integrals 
of motion strongly constraints the system dynamics and 
prevents it from reaching thermal equilibrium and explor- 
ing all possible states having the same average energy and 
particle number (grand canonical Gibbs ensemble). 

By using the reduced density matrices for this class 
of initial states, we can show that nonergodicity implies 
that the asymptotic correlation functions can be effec- 
tively described by an ensemble that assigns a mode- 
dependent temperature to each eigenmode. This is pre- 
cisely the physical content of the generalized Gibbs en- 
semble (GGE) [JJ. We have illustrated our method by 
analyzing the the quantum Ising and the XX spin chain 
models, both of which exhibit fermionic quasi-particles. 
In Sec. |III[ the Luttinger model, which exhibits bosonic 
quasi-particles, was analyzed by the same method. 

For initial states lacking the lattice translational sym- 
metry, we have shown the connection between nonergod- 
icity and the fact that eigenmode correlations in the ini- 
tial state become dominated by diagonal correlations (i.e. 
quasi-particle occupations) as the thermodynamic limit 
is approached. By direct numerical calculation and phys- 
ical reasoning, we have argued that this property should 
hold true for quantum quenches in which a physically 
sensible potential (e.g. a trap) that scatters the quasi- 
particles is suddenly removed at t = 0. However, at 
present we are unable to provide a mathematically rigor- 
ous proof of this fact (cf. Eq. 43), although in all cases 
that we have examined so far, it appears to hold true. 

Furthermore, using method reported here, we have 
been able to analytically shed light, for a much broader 
class of exactly solvable models and initial states than 
considered so far [SJ [H] , on the conditions under which 
the generalized Gibbs ensemble is expected to apply. 
Thus, our results extend the validity of the GGE con- 
jecture to a much broader class of quantum quenches. 
Our method also explains the special role played by the 
quasi-particle occupation operators as the set of integrals 
of motion required to construct the GGE. The noner- 
godic behavior of the correlation functions found here is 
entirely explained by the dependence on the expectation 
value of such integrals of motion alone. 

In future studies (26] [27] , it will interesting to under- 
stand how these results relate to the generalized eigen- 
state thermalization hypothesis discussed in Ref. [15] . 
We will also apply our methods to understand the con- 
ditions under which the asymptotic state can become ar- 
bitrarily close to a thermal state [2B]. The latter study 
unveils further interesting connections between the GGE 
and quantum Information theory [261 . 
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Appendix A: Eigenmodes of the quantum Ising and 
XX chains 



like crj in the quantum Ising and XX chain models de- 
pends on the applicability of Wick's theorem in the 
t — > +00 limit to multi-point correlation functions. Thus, 
we shall first tackle this problem by time-averaging the 
correlation functions of finite systems prior to taking the 
thermodynamic limit. Let 



Cg\xi,xj,t) = (O(x,t)O(0,t)), 



(Bl) 



for local operator like 0(x) = ^ fe <p k (x)f k . Its time av- 
erage is defined as: 



The Hamiltonian of the XX and Quantum Ising chains 
introduced in the main text can be brought to diagonal 
form by means of the Jordan- Wigner transformation: 

al = l-2fjf h a; \[:\ 2/; /,:/,. (Al) 

j<i 



C%> ( Xi , Xj ,t)= lim 1 / dt C% ] (x { , Xj , t) . (B2) 



(A2) 



with {.fi, fj} = Sij, anti-commuting otherwise. 

For the Quantum Ising chain, assuming periodic 
boundary conditions, that is, 



fj 



1/2 



(A3) 



with xj = j and 
f{k) =cos(^(fe)/2)7(fc) + isin(^(fc)/2) 7 t(-fe), (A4) 

where tan# g (fc) = sinfc/(cos/c — g). These two transfor- 
mations render the Hamiltonian of the Quantum Ising 
chain, Eq. [2j diagonal: 



7 t(fc) 7 (fc)-- 



(A5) 



where e g (k) = U\J 1 + g 2 — 2g 2 cos k. 

For the XX chain we shall assume an open ended chain 
where 



f> 



1/2 



^sin/cx-, f(k), 



(A6) 



L + 

with k — jizj, j,rn= 1, . . . , L, yields: 

H = Y J <k)f\k)f{k), (A7) 

k 

where e(k) = —Jcosk — h. 

Appendix B: Time averages and Wick's theorem 

As mentioned in the main text, the calculation of 
asymptotic correlation functions of non-local operators 



(2) 

A priori, the time average of Cq [Xi , Xj , t) followed by 
the thermodynamic limit yields 



D { o'( Xi , Xj ) =^2(p* k (xi)(p k (xj)N (k), 



lim C%\xi,Xj,t), 

t — >- + oo 



(B3) 



where the last limit is taken after the thermodynamic 
limit. However, taking the thermodynamic limit after 
time averaging is a subtle procedure. For instance, for 
the four-point correlation, 

^ O (xi, Xj , x m: x n ,t^j {C^ (x^ , £)0"^ (xj , i) 



x 0(x m ,t)0(x n ,t)), 



it yields: 



(Xi : Xj : X m , X n: t^j ^ ^ A kk > (^Xj , Xj , x m , x n ) 



k,k' 



x [N (k)N (k') - |G (M')| 2 ] , (B4) 



where 



A k ^(xi,Xj,x m ,x n ) = ip* k (x t )ip k (x n )(p* k ,( X: j)ip k ,(x m ) 

- tpl(xi)(p k (x m )(pl,(xj)(pk>{x n ). (B5) 

Since we have assumed that the square of the functions 
<fk(xi) is normalized to system size, i.e. \ip k (xi)\ ~ 
0(L -1 / 2 ), it follows that Ak,k'{xi,Xj,x m ,x n ) ~ 0(L~ 2 ), 
which is required to obtain a finite result in the L — > +00 
limit given the presence of the double sum over k and k' . 
However, we note that 



(.Xi? Xj , X mi X n ,t^) ^ Cq {Xi,X ni t^J Cq (x^,x m ,t) 



— Cq ^(Xi, X m , t) Cq ^(Xj, x n , t). (B6) 

and thus, if we also identify D^(xi,Xj,x m ,x n ) — 

lim t _j. +00 Cq (xi, Xj, x m , x n , t) with its time average, 
Wick's theorem appears to be violated for t —} +00 in the 
thermodynamic limit as D^(xi,Xj,x m ,x n ) will a priori 



11 



depend on all quantum correlations between the eigen- 
modes described by Go(k,k') (cf. Eq. B4|. However, 



D^ 2 \xi, Xj), which we identified with C( 2 )(ie.;, Xj, t), only 
depends on No(k) = Go(k,k). This has implications for 
the calculation of correlation functions of non-local op- 
erators in the t —> +00 limit. Nevertheless, as it was 
discussed in the main text, the eigenmode correlations, 

|G (fc, k')\ 2 -> [N Q {k)f 8 k ,„ + AR (k, k'), (B7) 

as the thermodynamic limit is approached. Therefore, 
the term involving these correlations in Eq. (B4), be- 



comes approximately equal (after neglecting AR (k,k') 
to 

-^2 A k,k(xi,Xj,x m ,x n )NQ(k), (B8) 



which is manifestly of 0(L 1 ) as L — >• +00. On the other 
hand, the term 



A k>k >(xi,Xj,x m ,x n )Na(k)No(k') 



(B9) 



is of O(L ) as L —> +00. However, Eq. (B9) equals the 



anti-symmetrized product of the time averaged two-point 
correlation functions, C( 2 ~>(Xi,Xj). Thus, the fact that 
the contribution of non-diagonal correlations becomes 
negligible in the thermodynamic limit justifies the proce- 
dure to taking the time-average followed by the thermo- 
dynamic limit. 

We can try to extend the above result to higher or- 
der (i.e. 3-point, etc.) correlation correlation functions. 
However, it is more convenient to take a shortcut. As 
discussed in the previous paragraph, we can identify the 
time average of two-point correlation functions of local 
operators like 0(x) with the t — > +00 limit of the same 
correlation function in the thermodynamic limit. These 
two point correlation function are the building blocks for 
computing with multi-point correlation functions or cor- 
relation functions of non-local operators like a x as we can 
apply Wick's theorem first and then let t —> +00 using 



M 



lim c 1 (t)c 2 (t)- 

t— >+oo 



•CmW = JJc 4 (+oo) 



(BIO) 



where Ci(+oo) = lim t _j. +00 a (t) since the t — > +00 limit 
of two-point every correlation function (ci(t) in the ex- 
pression above) exist and it is given by the its time aver- 
age followed by the thermodynamic limit. 



Appendix C: Non-local correlations in the quantum 
Ising and XX chain 



for t +00 can be computed first using Wick's the- 
orem and letting t — > +00 in the resulting expres- 
sion. To this end, it is convenient to write af(t) — 

Mt) = f\{t) + m) = (jAr[) E sin ^ [ eHk)t/h 



x p{k) + e~ ie{k)t/h f(k)], (C2) 
Bi(t) = f}(t) - Mt) = (^ — J J2 sinkx i [ e * 



Je(k)t/H 



x /t(fc)_ e -"W*/ft f(k)]. 



(C3) 



where the mode expansions are given for the XX chain 
(cf. e.g. Ref. [20] for the corresponding expressions for 
the quantum Ising chain). Hence, 



Cg>(x i ,x j ,t) = (B i (t) 



n Mm® 

i<i<j 



Mt)), (C4) 



where we have assumed that Xi < xj without loss of gen- 
erality. To evaluate the expression above we use Wick's 
theorem and take the limit t — > +00 of the resulting 
expression only after taking the thermodynamic limit. 
Moreover, since 



(Mt)A 3 (t)) ->5ij, 
(Bi(t)Bi(t)> -5ij, 



the correlation function 



■^XX } ) lim CxX (*^Z j *^ J ') 

t— 



can be written as a Tocplitz determinant 



(C5) 
(C6) 



(C7) 



ao ai • • • a_ n+ i 
Oil ao ■ ■ ■ a_ n+ 2 



Q-n-1 dn-2 



«0 



(C8) 



(2) 

where a,i-j+i = — D^g{Xi, Xj), and 
D% ) B (x i ,x j )= lim (Ai(t)Bj(t)) 



t— f + OO 



2 J £<p* k (x i )v k (x j )N (k), (C9) 



For the XX chain, the momentum distribution of the 
hardcore bosons can be obtained from the expression: 

1 \ - 



S(q,t) 



2L 



5 ^-*i)cW(i il x Jl t) 1 (CI) 



where <p k (xi) = \/2/(L + 1) sinkxi and the thermody- 
namic limit L — > +00 at finite lattice filling, hq = N/L, 
is implicitly understood. Note that, in this limit, the 
actual boundary conditions (open or otherwise) are irrel- 
evant. 
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